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SUMMARY 


This  research  presents  two  initialization  methods  in  the  form  of  optimal  filter 
equations  suitable  for  adiabatic  and  diabatic  numerical  weather  prediction  and  also  an 
optimal  system  of  equations  for  numerical  weather  prediction  involving  four  prognostic 
equations,  the  optional  mixing-ratio  continuity  equation,  and  one  diagnostic  equation. 

The  first  initialization  method  commensurable  with  the  hydrostatic  assumption  is 
a  practical  generalization  of  Fjortoft’s  results  and  docs  not  suffer  from  some  weakness¬ 
es  of  the  balance  cquation-particularly  in  the  anticyclonic  case  analyzed  by  Elsacsscr. 
It  also  significantly  improves  the  effective  prediction  time  scale  and  can  readily  replace 
computer  programmed  models  based  on  the  balance  equation  with  a  resulting  improve¬ 
ment  of  global  weather  prediction. 

The  second  initialization  method  consisting  of  three  filter  equations  in  the 
(x,  y,  z,  t)-system  does  not  presuppose  hydrostatic  equilibrium  and  does  not  encounter 
the  problem  of  hyperbolicity  and  is,  therefore,  applicable  to  smaller  scales.  The  third 
filter  equation  may  simultaneously  replace  the  hydrostatic  one  within  an  optimal  prog¬ 
nostic  system  containing  the  horizontal  equations  of  motion,  the  continuity,  and  the 
thermodynamic  equation  in  their  invariant  forms.  The  relaxation  procedure  for  solving 
the  nonhydrostatic  filter  differential  equations  is  no  more  difficult  than  the  numerical 
solution  of  the  general  balance  equation  in  conjunction  with  the  co-equation. 

Finally,  the  study  throws  new  light  on  problems  such  as  ergodieity,  long-range 
predictability,  gridscale  diffusion,  surface  friction,  and  boundary  conditions. 


FOREWORD 


This  research  was  started  in  1965  and  is  a  persona,  effort  except  for  editing,  typ¬ 
ing,  and  printing.  Some  of  its  results  were  already  available  in  1966.  In  the  meantime, 
the  author  tried  to  perfect  the  analysis  and  to  put  it  in  the  right  perspective.  It  is  of 
significance  that  the  filter  methods  developed  for  the  hydrostatic  and  les«  restricted 
(x,y,z,t)-system  associated  with  a  new  set  of  prognostic-diagnostic  equations  are  not 
only  new  and  optimal  but  have  a  useful  purpose  and  include  opcrativencss  as  well.  The 
author  would  like  to  thank  Colonel  John  R.  Oswalt,  Jr.,  who  endorsed  the  preparation 
and  preliminary  publication  of  this  work  in  the  form  of  a  comprehensive  USAETL  Re¬ 
search  Note. 
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THE  DERIVATION  AND  POTENTIAL  ,)F  NEW  FILTER 
EQUATIONS  FOR  NUMERICAL  WEATHER  PREDICTION* 


1.  introduction.  The  Global  Atmospheric  Research  Program  (GARP)  spoil- 
"  'red  by  the  International  Council  of  Scientific  Unions  and  the  World  Meteorological 
^  ganWation  is  aimed  during  the  next  decade  at  gaining  additional  insight  about  the 
•  '  of  the  atmosphere  so  that  !  becomes  feasible  to  predict  the  weather 

;i '  ide  for  periods  r*‘  at  least  2  weeks. 

The  first  requirement  to  meet  this  goal  is  the  establishment  of  a  worldwide 
observation  system  involving,  advanced  satellites  providing  for  interrogating,  recording, 
and  locating  on  a  realtime  basis  and  for  measuring  of  cloud  cover.  The  satellites  will 
carry  passive  microwave  radiometer,  occultation  and  refraction  sensors  for  the 
determination  of  vertical  profiles  of  temperature,  water  vapor,  and  density  from  which 
constant  pressure  heights  could  be  determined.1 

The  second  requirement  is  to  acquire  a  better  understanding  of  some  physical 
processes,  such  as  the  flow  of  heat  and  moisture  near  the  tropical  waters,  which  have 
to  be  incorporated  in  a  long-range  atmospheric  prediction  scheme.  The  processes  of 
air-sea  interaction  and  related  research  efforts  are  the  subject  of  the  Barbados 
Oceanographic  and  Meteorological  Experiment  (Bomex)2  which  is  the  first  experiment 
of  the  United  States  portion  of  GARP. 


The  third  requirement  is  the  development  of  computers  capable  of  processing  the 
tremendous  amount  of  data  and  simulating  the  extremely  complicated  atmospheric 
generation  process.  The  Bomex  repi  "t  states  that  a  computer  now  being  designed  is 
expected  to  have  100  io  1000  times  the  capacity  of  present  models. 

Finally,  more  sophisticated  prccnction  models,  including  nonadiabatic  and 
'fictional  effects,  some  of  these  at  least  in  a  parametric  form,  will  have  to  be 


developed.  As  Charncy  stated  in  195 13  and  since  proved  by  others,  condensation 
phenomena  are  the  simplest  to  introduce  although  much  remains  to  be  done  in  this 
respect.  The  same  is  true  for  long-wave  radiational  effects  which  are  prognostieally 
interrelated  with  humidity  forecasts.  Very  difficult  to  incorporate  is,  of  course,  the 
turbulent  transfer  of  momentum  and  heat. 

Within  the  context  of  BOMEX,  a  “Hemispheric  Model  Study”  is  to  be  conducted 
by  Pandolfo,  while  Charncy  acts  as  principal  investigator  of  a  “Theory  of  Large-Scale 
Atmospheric  and  Oceanic  Processes.” 

As  to  the  utility  of  improved,  short-range  predictions  including  those  of 
precipitation  and  of  useful  long-range  forecasting,  it  should  be  mentioned  that  the 
National  Academy  of  Sciences  —  National  Research  Council4  estimated  the  potential 
savings  resulting  therefrom  as  approaching  2Vi  billion  dollars  annually  for  the  United 
States  alone. 

As  Thompson5,  6  and  Novikov7  have  pointed  out,  and  recognized  again  in  the 
panel  discussion  on  atmospheric  predictability  during  the  Golden  Anniversary  Meeting 
of  the  American  Geophysical  Union  in  April  1969,  and  also  evident  in  Sections  4 
through  7,  herein,  the  initial  state  of  the  atmosphere  is  of  considerable  importance  for 
long-range  and  sometimes  even  for  short-range  predictions  in  cases  of  rapid 
development.  F.G.  Shuman  states:  “We  are  aware  of  many  weaknesses  in  our 
initialization  procedures,  principally  the  lack  of  the  full  meteorological  balance 
between  the  windc  and  pressure  field.  This  problem,  however,  is  perhaps  even  more 
difficult  than  those  encountered  in  developing  the  prediction  model.  The  balance 
problem  should  properly  be  approached  as  part  of  the  analysis  problem.  At  NMC 
(National  Meteorological  Center)  we  are  making  a  start  on  th '  problem,  but  it  will  be 
some  time  before  we  have  a  general  solution  to  it.  We  expect  to  gain  a  great  deal  from 
other’s  work  on  the  problem.”8 

After  i  critique  of  present  initialization  schemes,  this  report  presents  more 

powerful  initialization  methods  in  the  form  of  optimal  filter  equations  for  numerical 

weather  prediction  which  do  not  suffer  from  some  essential  weaknesses  of  the 

3J.G.  Chamey:  “Dynamic  Forecasting  by  Numerical  Process,”  Compendium  of  Meteorology,  American 
Meteorological  Society,  Boston,  Massachusetts,  1951. 

4National  Academy  of  Sciences  -  National  Research  Council:  “Useful  Applications  of  Earth-Oriented  Satellites," 
Meteorology,  National  Academy  of  Sciences,  Washington,  D.C.  1969. 

5P.D.  Thompson:  Numerical  Weather  Analysis  and  Prediction.  The  MacMillan  Company,  New  fork,  1961 . 

6P.D.  Thompson:  Uncertainty  of  Initial  State  at  a  Factor  in  the  Predictability  of  Large-Scale  Atmospheric  Flo te 
Patterns,  Telius,  1957. 

7E,A.  Novikov:  “Toward  a  Formulation  of  the  Problem  of  the  Predictability  of  Synoptic  Processes,”  U.S.  Joint 
Publications  Research  Service,  New  York,  August  1960. 

*F,G.  Shuman:  “A  Multi  level  Primitive  Equation  Model,”  U.S.  Department  of  Commerce,  Environmental 
Sciences  Service*  Administration,  Weather  Bureau,  National  Meteorological  Center,  November  1965. 
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I  so-called  balance  equation,  i.c.,  they  are  more  accurate,  applicable  to  relatively  smaller 

f  scales,  and  suitable  to  include  convection  in  case  of  the  general  nonhydrostatic  filter 

process.  The  new  filter  equations,  one  of  which  simultaneously  provides  a  diagnostic 
'  equation  for  an  ultimate  prognostic  system,  can  also  be  profitably  used  for 

re-initialization  and  ire  fundamental  pertaining  to  the  inclusion  of  nonadiabatic 
:  processes  through  the  use  of  statistical  and  parametric  procedures.  The  general 

nonhydrostatic  filter  equations  also  offer  an  interesting  and  more  satisfying  approach 
as  far  as  the  use  of  numerical  relaxation  is  concerned.  Furthermore,  the  analysis 
throws  some  new  light  on  the  significance  of  initialization  for  long-range  pre¬ 
dictability,  aid  the  use  of  frictional  terms  and  boundary  conditions.  Finally,  the  new 
equations  offer  advantages  with  reference  to  the  computation  of  winds  from 
worldwide  data  obtained  through  the  use  of  advanced  satellites. 


2.  Thompson’s  Generalized  Filter  Approximations,  the  Balance  Equation,  and 
Related  I'  ion  Methods.  Thompson9  showed  that  it  is  necessary  and  sufficient  to 
filter  gravity  waves  out  if  the  total  time  derivative  of  the  divergence  is  omitted  from 
the  “divergence”  equation  obtained  by  applying  the  horizontal  divergence  operator  on 
the  equations  of  motion.  This  results  in  *he  diagnostic  relationship 


Vu  •  4—  +  Vv  •  -  k  •  V  x  fV  +  A*  'F  =  0  (1) 

3x  by  ' 

-> 

where  V  is  the  horizontal  nabla  operator;  V,  the  horizontal  velocity  vector;  k  a 
vector  directed  vertically  upward;  f  ,  the  Coriolis  parameter;  A1 ,  the  2-dimensional 
Laplace  operator;  and  'F  ,  the  isentropic  stream  function  in  the  (x,  y,  0 ;  t)  -  system 
where  0  is  the  potential  temperature.  Equation  (1)  is  valid  for  adiabatic,  nonviscous 
flow  and,  of  course,  under  the  hydrostatic  assumption. 

Omission  of  with  D  as  divergence  in  eq.  (1),  allows  the  establishment  of 

the  diagnostic  w -equation: 


=  £  Kvlv)  -V!(H 

+  •  *p)  • 


Here,  it  is  w  =  ,  p  the  density,  and  p  the  pressure,  while  K,  R,  and  B  arc  specific  func- 

dt 

tions  free  of  time  derivatives.  In  order  to  obtain  eq.  (2),  use  has  to  be  made  of  the  vorti- 
city  equation  and  the  continuity  equation  in  isentropic  coordinates. 

9P.  D.  Thompson,  “Numerical  Weather  Anilyeu _ ” 


Since,  according  to  Helmholtz'  theorem,  the  horizontal  vector,  V,  may  be 
expressed  as  the  sum  of  an  irrotationnl  and  a  nondivergent  vector  so  that 

V  =  Vx  +  k  x  7*  , 

curl  V  -  k  A2 1//  and  div  V  =  A1  x  =  +  |~  •  7p  j  . 


Under  hydrostatically  stable  conditions,  eq.  (1)  and  (2)  may  thus  be  solved 
for  ^  and  x-  For  prognostic  purposes,  either  the  vorticity  equation  together  with 
eq.  (1)  or  (2)  or  the  primitive  equations  together  with  the  continuity  equation 


A  in  =  _  111  (iiv  V 
(It  W  do  <l,v  v 


might  be  used.  In  the  latter  case,  instabilities  might  occur  under  certain  conditions. 


A  simplification  of  eq.  (1)  as  well  as  the  starting  equation  for  an  iterative  process 
involving  eq.  (I)  and  (2)  is  the  so-called  balance  equation 

f A2  \p  +  2  (\pxx  \pyy  -  \p;y  )  +  V\p  •  7f  =  A2  vl>  (3) 


which,  according  to  Rellich’s  theorem10  and  as  demonstrated  by  Bolin}1  is  of  the 
elliptic  type  if 


A2'!' 

f 


-jrvr  .  v*  > 


£ 

•> 


(4) 


This  condition  is,  according  to  Bolin}1  also  necessary  to  solve  eq.  (2)  as  an  elliptic 
partial  differential,  equation  under  the  assumption  of  a  hydrostatically  stable 
atmosphere. 

Apart  from  the  fact  that  eq.  (4)  may  not  be  fulfilled  in  the  case  of  strong 
anticyclones,  Bolin11,12  has  pointed  out  that  the  most  serious  approximation  in  eq.  (1) 

is  neglecting  the  time-dependent  term  div  V  through  which  all  gravity-inertia 

oscillations  are  eliminated  and  that  the  filtering  approximation  is  only  valid 
for  div  V  <  f  =  A2  ip .  As  mentioned  by  Phillips,13  fast  moving  and  nongcoslrophic 
gravity  waves  have  large,  horizontal  divergences  (div  V  ~  £ ),  and,  during  convection, 

10  F.  llcllich:  “Zur  ersten  KandwcrUtifgabc  l>rl  Monge-Ampereiwhfn  nifferentlalgleielningpn  votn  rllipllsrhrn  Typas; 
differentialgeomctrischc  Anwendungcn,”iWiiil6.  Ann.,  Vol.  107, 1933, 

11  II.  Itolin:  "Numerical  Forecasting  with  the  Barotroplc  Model, 'TeMut  VII,  1935, 1. 

12  It.  Itolin:  “An  Improved  Barotroplc  Model  and  Some  Aspect*  of  Using  (hr  Balance  Equation  for  Thrcc-dimeWsional 
Flow,"  Telhu  VIII.  1956, 1. 

13  N.  A.  Phillips:  Numnrtcnl  Weather  Prediction,  Advance*  in  Computer  Control  (or  Technique*),  Vol.  1,  Academic 
Press,  New  York,  1960. 
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div  V  >  f  . 


Although  Grossman14  anti  other  authors  attributed  a  considerable  improvement 
in  numerical  forecasts  lo  the  use  of  the  balance  cq,  (3),  Elsaesser1*  has  come  to  the 
conclusion  (hat  its  linearized  version  appears  to  be  the  optimum,  particularly  in  view 
ol  recovery  ol  the  pressure  height  (if  is  replaced  by  0  in  the  (x,  y,  p)  —  version 
ol  et|,  (3)),  less  computation  time  and  the  necessity  of  elliptizing  in  case  eq.  (4)  does 
not  hold.  Elsaesser  explains  the  differences  with  respect  to  previous  investigations. to 
the  near  hemispheric  scope  of  his  study  and  the  inclusion  of  large  areas  of  strong 
anticyclonie  shear  equatorward  of  the  polar  jet  stream. 

Ilinkelmann16  has  demonstrated  that  the  stream  function  obtained  through  eq. 
(3)  fails  during  strong  anticyclonie  activity,  and  a  belter  result  would  be  achieved  by 
applying  the  gcoslrophic  stream  function.  Monin17  has  emphasized  that  eq.  (3)  would 
be  equivalent  to  an  equation  in  which  the  quadratic  terms  would 
read  -  02xy )  as  *ar  as  l*lc  rtoa*<'  ^,<!ory  >s  concerned.  Arnason18  ,  who 

neglects  Ihcir  vertical  advcction  terms  in  the  equation  of  motion,  arrives  at  similar  and 
practically  equivalent  results. 

A  method  of  initialization  using  finite  difference  techniques  and  employing  the 
balance  condition  -j^-  div  V  =  0  has  been  developed  by  Miyakoda  and  Moyer”  .  It  is 

simpler  from  the  computational  point  of  vie;/  than  the  classical  method  of  solving  the 
balance  equation  together  with  the  (^-equation,  and  numerical  results  obtained 
through  the  use  of  a  simplified  model  agree  rather  well  except  for  a  slight  displacement 
of  the  vertical  velocity  pattern.  Moreover,  this  method  would  be  relatively  more 
advantageous  if  the  initial  vertical  velocity  bud  an  insignificant  influence  on  the 
meteorological  evolution.  This  is,  however,  only  true  in  case  of  weak  development  and 
also  due  lo  the  use  of  the  balance  equation  (as  will  be  shown  later).  Better  filter 
equations  are  indeed  necessary  for  forecasts  exceeding  5  days  after  which,  according  to 
the  experience  gained  by  Mint/,  and  Krishnamurti,  as  quoted  by  the  authors,  rather 
exact  initial  conditions  become  important. 


M  (».  I’.  Crewman:  "A  Study  of  Numerical  Forecasting  F.rrors,”  Monthly  Weather  Review,  July  1957. 

15  It.  W,  Kllsacsscr:  “Comparative  Tent  of  Wind  Law*  for  Numerical  Weather  Prediction,"  Monthly  Weather 
Review,  Vol.  92,  No.  5,  1 968. 

16  K.  Ilinkelmann:  “I teller  die  Finbrzichunff  dlvergentcr  Wlmlkomponenten  in  VorhersaRrmodclle,"  flrrirhtc  det 
deuttchtm  Welturdiemtei  Nr.  30  (Hand  5)  Had  Kiiwingun,  1957. 

17  A.  S.  Monln:  “Izcmcncnila  davlcniia  v  baroktlnnol  atmosfnrc  tzv.  Akad.  Nauk,"  ter.  fr'ofa.,  no.  4, 195R. 

18  (i.  A  meson:  “Illghcr-Ordcr  Ocoslrophic  Wind  Approximations,  " Monthly  Weather  Review,  Vol.  90,  No.  5, 
May  1962. 

19  K.  Miyakoda  and  R.  W.  Moyer:  “A  Method  of  Initialization  of  Dynamical  Wrathcr  ForeeasllnR,"  Geophysical 
Fluid  Dynamics  Laboratory,  environmental  Sciences  Services  Administration,  Washington, October  1966. 
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According  to  Phillips,30  noise,  through  an  initial  divergence  implied  by  the  usual 
geostrophic  theory,  is  greatly  suppressed  if  the  initial  data  is  used  for  the  primitive 
equations.  However,  the  suppression  of  high-frequency,  gravity-inertia  waves  will  only 
be  very  effective  if  the  motion  is  not  strongly  ageostrophic:  “The  extent  to  which  the 
method  will  prove  useful  in  low  latitudes,  or  in  high  latitudes  when  div  V  ap¬ 
proaches  i  in  magnitude,  can  perhaps  be  answered  only  by  experiment.”30 

A  method  (similar  to  the  method  of  Miyakoda  and  Moyer)  by  Nitta  and 
Hovermale31  involves  an  actual  iteration  of  forward  and  backward  forecasts  around 
the  initial  time  with  the  Euler-backward  time  difference  and  yields  acceptable 
rotational  but  unsatisfactory  divergent  wind  components.  According  to  Nitta  and 
Hovermale,  “the  question  remains  unanswered  whether  o1  not  the  technique  in  its 
present  design  is  accurate  enough  to  be  of  practical  value. . . .  Perhaps,  further 
diagnosis  . . .  will  lend  some  insight  in  this  direction.” 

This  diagnosis  is  readily  available  since  the  method  corresponds  to  the  filter 
conditions 


{&)  ■(#) 

'  't=o  %  'tn-o 


However,  this  filtering  leaves  the  factorial  term  in  the  equations  of  motion 

undetermined  so  that  in  addition  =  0  has  to  be  postulated  which  is  essentially 
equivalent  to  f-A  t~-\  =  0  and,  thus,  slightly  less  satisfactory  than  the  balance 

V  dp/to=0 

equation  together  with,  the  corresponding  u  -equation. 


The  co -equation  corresponding  to  the  balance  equation  has  been  investigated  by 
many  authors.  Krishnamurti33  uses  the  vorticity  equation  including  frictional  terms 
together  with  the  balance  equation  and  the  w  -equation  in  the  (x,  y,  p,  t)-system  for 
diagnostic  studies  of  weather  systems  with  a  term  for  latent  heat  included.  Vukovich 
and  Chow30  employ  the.  balance  u -equation  in  order  to  estimate  the  effect  of 
long-term  diabatic  heating.  Considerable  experience  in  this  regard  has  also  been  gained 
by  Smagorinsky  and  collaborators.  (Some  of  their  results  will  be  discussed  in  Section 
5.) 

30  N.  A.  Phailp*:  "On  the  Problem  of  Initial  Data  for  the  Primitive  Equation*,”  Teliut,  VoL  12,  May  1960. 

31 T.  Nitta,  and  J.  B.  Hovenndc:  “On  Anatyaia  and  Initialisation  for  the  Primitive  Focecaat  Equation*," 
Technical  Memorandum  WBTM  NMC-42,  National  Meteorological  Center,  Suitland,  Maryland,  October  1967. 

33  T.  N.  Kriahnamurti:  “Diagnostic  Studies  of  Weather  Systems  of  Low  and  High  Latitudes  (Romby  Number  1),” 
Department  of  Meteorology,  Unhrenitv  of  California,  L  oa  Ai^elea,  California.  30  Nov.  1964  -  30  Nov.  1966. 

°  F.  M.  Vukovich  and  Chi  F.  Chow;  "Research  in  Numerical,  Dynamical,  and  Operational  Weather  Fere  out 
Procedure*,"  Research  Triangle  Institute,  Research  Tri*i*ie  Park,  North  Carolina,  April  1968. 
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“  l£j3h-  z  ; 


In  conclusion^  it  should  be  emphasized  that  the  shortcomings  of 
equation,  i.e.,  vj  =0  and  A20  +  -y>  0  where  0 


the  balance 
denotes  the 


geopotential,  are  necessarily  reflected  in  the  more  explicitly  named  balance 
cj -equation.  Superior  filter  equations  which  do  not  suffer  under  these  restrictions  wili 
simultaneously  result  in  an  overall  improvement  of  the  initialization  process  and,  thus, 
in  better  short-  and  long-range  forecasts,  as  well  as  in  the  elimination  of  high-frequency 
meteorological  noise  and  in  stabilization. 


3.  Fjortoft’s  Filter  Equations.  The  application  of  the  filter  conditions 

(ft.  #1.  ■ 0 


(5) 


with  reference  to  the  equations  of  motion  in  the  (x,  y,  p,  t)-system  by  Fjortoft24  leads 

to  two  prognostic  filter  equations  with  the  local  time  derivatives  and  —  ~~  in 

01  ox  ot  dy 

which  0  designates  the  geopotential.  Together  with  the  adiabatic  thermodynamic 
equation  involving  and  the  continuity  equation  -~  +  +  -|~  =  0  where 

do  c)  a 

to  =  -g*-,  these  are  four  equations  for  the  unknowns  u, w,  and  —  .  As  the  ellip- 

ticity  criterion  under  which  a  solution  with  a  suitable  relaxation  factor  a  of  the 
same  sign  can  be  achieved,  Fjortoft  lists  the  conditions 

(f  **^**^-^*{»**f  *){*>,  -F0p)-(v0p)  >0  (6) 

where  F  = 


xy 


0, 


*p 


R 

h5<5 

and 

+  0 

Yxx 

0 

*xy 

0 

0 

vyp 

>  0. 

(7) 

0  -  F0 

p 

Approximately,  the  conditions  may  be  stated  as  f2  +  A2  0  >  0  which  means  a  considerable 

fj 

improvement  in  the  anticyclonic  case  when  compared  with  +  A2  0  >  0  applicable  in 


24  R.  Fjortoft:  "A  Numeric*!  Method  of  Solving  Cert* in  Partial  Different! *1  Equations  of  Second  Order,"  Fremlsgt 
1  Videnakapa-Mote,  February  1962. 
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connection  with  the  balance  equation,  and  $  -  V<pp  >  0  or  a  >  0  with  a  as  the  static 

stability. 

Because  Fjortoft  has  no»  developed  superior  diagnostic  filter  equations  based  on 
the  same  filter  conditions,  his  method  does  not  appear  to  have  been  practically 
utilized. 

As  to  the  application  of  a  relaxation  factor  of  changing  sign  in  the  hyperbolic 
case,  no  useful  results  may  be  obtained  since  the  simple  model  of  the  continuity 
equation  does  not,  or  not  sufficiently,  hold  and  should,  therefore,  not  be  adapted  to 
the  relatively  more  sophisticated  thermodynamic  equation.  Apart  from  the  fact  that 
no  accurate  results  can  be  expected,  particularly  in  the  case  of  pronounced 
hypcrbolicity,  the  convergence  of  the  relaxation  process  iias  been  extremely  slow 
according  to  Elsaesscr.25 

4.  Optimal  Filter  Equations  in  the  (x,  y,  p,  t)-System.  For  the  sake  of  simpli¬ 
city,  we  omit  in  the  following  paragraph  the  use  of  map  scale  factors  as  has  been  done 
by  Fjortoft.  Inclusion  of  these  factors  docs  not,  however,  present  any  fundamental 
difficulties.  Moist-adiabatic  and  frictional  terms  are  also  neglected  at  present.  Thus, 
we  proceed  from  the  following  system  of  equations  which  is  familiar  to  the  reader: 


du  dz  ,  . 

df  -g  3x  +  fv 

(8) 

dv  dz  , 

57 “  fu 

(9) 

|l  +  i"  =  0 

dx  dy  dp 

(10) 

dT  RT  _  n 

dt  Cpp 

(11) 

u  *  i  3z  .  1  _  RT  .  _  dJz  .  1  dz 

or  alternatively,  with  -r—  -  —  -  —  and  a  -  rr  +  — —  -r— 

dp  gp  gp  OF  *P  °P 


-jr— •  +  U  +  V-|-  ^  +  OGO  =  0  . 

dt  dp  dx  dp  dy  3p 

After  applying  the  filter  operations  of  ea.  (5).  we  arrive  at 

d  dz  d  ,, .  .  _  n 
®  dt  dx  3t  ^ 

bH.  W.  EIlMewer. 


(12) 
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attftfi iiasa  i’i hata 


K  nwwwn-.*.-,*  r  V  -n; .  ....  .  •  ..*«•  .  ir-nr*  ...  rf^pmfn'm. • .  ,  :  -r(,  ■jrm>";  1  <'■ 


d  9z  .  d  j&  y  « 

S¥ay  +  lt(fu>  ~  0 

and  more  explicitly  at 


9  dz  .  8  9 i  ,  8  3z  |  3  3z 
3t  8x  U  3x  3x'  V3y’3x  w  3p  3x 


v/  3f  ,  3f\ 

"g(U3x  v3yj 


l(g|-fa)  =  o 


3  3z  .  3  3z  3  3z  |  3 

3t  3'y  11  3x  8y  v  3y  3y  u  3p 


3z  .  u  /  3f  3f  \ 
3y  g\u3x  v8yj 


+t(-*S + *)  ■  »• 


(13) 


(14) 


In  order  to  make  the  two  filter  equations  (13)  and  (14)  consistent,  we  eliminate  tlie 
time  derivatives  by  applying  the  operator  ~  on  eq.  (13)  and  —  on  eq.  14  and 
subtracting.  This  yields,  after  omitting  irrelevant  terms,  the  diagnostic  filter  equation 


f2  /3v  3u\  3u  blz  ,  ( 

g  \3x  “  by)  -  by  W  V 


3u  3v\  b2z  3v  32z 
3x  ~  by)  bxby  bx  dy5" 


3co  3lz  3 co  3Jz 


3x  3y3p  3y  3xdp 


f  A3  z. 


After  we  differentiate  eq.  (13)  and  eq.  (14)  pertaining  to  x  and  y  respectively  and  add 
the  results,  we  obtain  as  a  corresponding  dynamic  filter  equation 


-A2z  =(—  +  asz^  —  il*_  +  3w  + 

t  \ g  /  3p  \3x  3x3p  3y  3y3p, 


A  3Jz  .  3u  32z  , 
'jbxby  +  37  37*  + 


(15) 


+  il/l* .  2f , 


(16) 


(19) 


f*  fi 

or  approximately  —  +  A2  2  >  0  which  is  considerably  less  restrictive  than  —  +  — tj-  >  0 
g  g 


of  the  balance  eq.  (1). 


The  ellipticity  condition  for  the  new  w  -equation  (17)  is  that  the  quadratic  form 


o 

0 


0 


0 


1  3*2 

2  dxdp 

_L  d22 
2  dydp 


IJl*.  *  Jll  i!+4»2 

2  dxdp  2  dydp  g  . 


(20) 


be  positive  definite  which  results  in  the  inequalities 


a  >  0  and  a 


>  0  . 


(21) 


In  order  to  insure  ellipticity  in  the  cyclonic  case  if  a  >  0  and  the  inequality  (21) 
does  not  hold  a  priori,  a  has  to  be  adjusted  to  a  slightly  positive  value  oc  so  that 
the  relation 


+  A 


(w)’] 


=  e,  >  0 


(22) 


with  e,  as  a  small  quantity,  is  fulfilled  which  is  essentially  in  agreement  with  proce¬ 
dures  applied  by  Smagorinsky26  and  Krishnamurti.27  In  anticyclonic  conditions, 


P_ 

g 


+  A2  2  <  0  ,  the  height  field  has  to  be  smoothed  so  that 


P 

g 


+ 


AX  =  6j  >  0 


(23) 


with  e2  as  a  small  quantity.  The  smoothing  results  simultaneously  in  more  favorable 


*J.  Smagoriniky:  “On  the  lnd"**o*i  of  Mokt  Adlib* tie  Processes  in  Numerical  Prediction  Model*,'*  Btrkhusdet 
DiuUchtA  WttttrtlientUi  Nr.  38  (t5ind  S)  Bad  KMnfen,  1957. 

®T.  N.  Krirfmamurti. 
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values  o.  ,  »  an^  a^T'  so  that 

A  ’  3x3p  3y3p 

-ff  *  *'-)-i|(&,)'  *  (-&)'! 


>  e3  «  0 


reasonably  holds.  Alternatively,  in  the  anticyclonic  case,  terms  involving  products  of 
second-order  derivatives  in  eq.  (17)  and  (18)  are  to  be  neglected  in  a  numerical 
integration  by  over-relaxation  for  all  points  for  which  the  inequalities  (21)  and  (19), 
respectively,  are  not  valid.  It  should  be  pointed  out  again  that  the  inequalities  (22) 
through  (24)  are  necessary  in  case  of  initial  hyperbolicity  in  order  to  adapt  w  in  eq. 
(12)  to  be  the  greater  scale  of  the  horizontal  divergence  of  the  filter  equation  (16) 
which  is  rather  a  smoothed  .  Hyperbolicity  can  only  be  eliminated  or  is  not 

present,  respectively,  if  the  hydrostatic  filter  condition  is  not  utilized  for  the 
derivation  of  a  continuity  equation  such  as  (10),  i.e.,  if  w  and  ijfi.  are 
maintained. 

Lateral  boundary  conditions  for  eq.  (18)  are  =  -f-\  both  for  a  physical 

boundary  and  that  of  a  hemispheric  grid.  For  the  u -equation,  approximate  lateral 
boundary  values  «L  may  be  determined  by  applying  the  regression  eq. 
AawL  =  -  X(x,y,p)wL  toeq.  ( 17)  so  that  an  ordinary  differential  equation  results 
(mentioned  by  Reuter1®  as  well  as  by  Gliassen  and  Klcinschmidt19  and  generally 
preferable  to  Charney’s  advective  model30).  In  this  respect,  eq.  (17)  may  be  further 
simplified  to  be  consistent  with  geostrophic  approximations.  Ad  usual,  o>0  =  0  at 
the  top  of  the  atmosphere,  and  wp  »  -  p  V  •  Vh  <f>,  with  </>,  (x,y)asthe  geopo¬ 
tential  of  the  ground.  ° 


It  is  characteristic  of  the  (x,  y,  p,  t)-system  that,  initally,  u  and  v  are 
required  only  at  the  lower  boundary  for  the  computation  of  <up  and  that 

Uj  and  v2  can  be  easily  determined  from  the  Poisson  eq.  A1  x  =  -  » to  be 

solved  at  each  p-level.  Due  to  the  present  insufficiency  of  observations  and  restriction 
to  relatively  few  levels,  the  emphasis  has  been  essentially  on  a  rather  accurate 
solution  eq,  (J8)  which  may  be  followed  by  a  determination  of  a)*1*  from  eq 

(17)  by  numerical  relaxation. 


*  H.  Reuter  “Zum  gegenwaertigen  Stand  der  numcriKhen  Wcttervorhrroagc,”  Oitemichischei  Ingenieur-Archiv, 
Band  X,  Heft  2-3,  Springer- Yeriage  Wien,  1956, 

®  A,  Eliaaaen  and  E.  Kieinachmidt:  “Dynamic  Meteorology, "  Handbuch  der  Phytik,  Bd.  XLV1H,  Springer  Verlag, 
Berlin,  Goettingen.  Heidelberg,  1957. 

XJ.  G.  Chamey. 


Alternatively,  4>  and  to  may  be  determined  in  conjunction,  and  the  problem 
of  finding  a  suitable  over-relaxation  coefficient  along  the  lines  of  Stuart  and  O’Neill31 
and  O’Brien”  would  practically  be  confined  to  the  w -equation.  This  numerical 
process  is  somewhat  complicated  by  the  necessity  of  finding  the  divergent  wind 
components  intermediately.  Fortunately,  the  J/-  and  co-equations  are  both  essentially 
linear  so  that  initial  guess  and  convergence  are  not  critical. 


The  inclusion  of  linear  frictional  terms  in  the  filter  equations  does  not  present 
difficulties.  Nonlinear  terms  would,  however,  cause  computational  instability  and 
would  also  be  incompatible  with  the  requirement  of  strong  smoothing  in  case  of 
convection  since  the  isobaric  divergence  equation  cannot  cope  with  situations 
involving  effective  static  stabilities  H2  <  0.  In  these  situations,  a  convective  tempera¬ 
ture  adjustment  such  as  that  developed  by  S.  Manabe33  has  to  be  made.  All  these  diffi¬ 
culties  and/or  adjustments  can  be  avoided  by  a  general  nonhydrostatic  system  of  filter 
equations  developed  in  Sections  6  and  7. 


5.  Ramifications  of  Optimal  Hydrostatic  Filter  Equations.  The  advantages  of 
cqs.  (15)  and  t\ 7)  in  relation  to  Thompson’s  equations  (1)  and  (2)  are  their  greater 
versatility  in  case  of  strong  anticyclonic  movements,  a  better  determination  of  winds 
including  the  elimination  of  the  weakness  of  the  balance  equation  with  respect  to 
anticyclonic  winds  ascertained  by  Elsaesser,34  and  the  agreement  with  Hollman’s 
analysis35  regarding  the  failure  of  the  balance  equation  in  strong  anticyclonic 
situations  and  the  use  of  the  geostrophic  instead  of  the  gradient  wind.  Fjortoft’s  filter 
conditions  (5),  as  well  as  ours,  involve,  of  course,  the  second  total  time  derivative  so 
that  standing  gravity -inertia  waves  of  the  form  ?"=A  sin  kxe*"11  sin  p  t  and  f"  = 


A  sin  kxsin  mz  sin  vt  are  preserved  and 


¥=  0.  Scale  inconsistencies  and  ad¬ 


justment  difficulties  inherent  in  Fjortoft’s  method  are,  however,  avoided  through  the 
establishment  of  diagnostic  equations  and  the  corresponding  imposition  of  Helmholtz’ 
theorem. 


Apart  from  a  better  determination  of  the  stream  function  4>  in  the  cyclonic 
case  as  well,  the  w -equation  (17)  also  permits  the  computation  of  stronger 
divergences/ convergences  and  vertical  velocities  because  z  >  A2  so  that  the 

vertical  influence  increases  relatively  to  the  horizontal  influence  in  a  solution  involving 

D.  #.  Stuart  tnd  T.  H.  R.  O’Neill:  “The  Over-Relaxation  Factor  in  the  Numerical  Solution  of  the  Omega 
Fetation,  Monthly  Weather  Review,  Vol  95,  No.  5,  May  1%7. 

aJ.  J.  O’Brien:  “The  Over-Relaxation  Factor  in  the  Numerical  Solution  of  the  Omega  Equation,”  Monthly 
Wither  Review,  Vol  96,  No.  2,  February  1968. 

33  S.  Manabe:  “Simulated  Climatology  of  a  General  Circulation  Model  vith  a  Hydrologic  Cyde,’’  Monthly  Weather 
Review,  Vol  93,  No.  12,  December  1965. 

*H.W.Ek_eaaer. 

*  G.  Hollmann;  “tlcbcr  die  Grenzen  der  geostrophiachen  Approximation  und  die  Einfuchrung  nichtgeoatroohiachcr 
Methoden,”  Beriehte  Jet  Deuteehen  Wetterdlenttet,  Nr.  38  (Band  5\  Bid  K iwngen,  1957. 
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a  (Ireen’s  function.  A  practical  demonstration  of  this  has  been  given  by  Smagorinsky36 

who  utilized  the  geostrophic  vorticity  equation  including  the  full  term 

and,  thus,  subconsciously  an  expression  approaching  equation  (16),  with  resulting  ver¬ 
tical  velocities  on  the  order  of  30  cm  sec'1 . 


The  importance  of  the  new  filter  equations  becomes  particularly  evident  if  we 
develop  the  absolute  vorticity  in  iscntropic  coordinates  in  a  Taylor  series  which  yields, 
in  view  of 


+f)  =  -(t»  +fKv 


(25) 


-  div  V, 
*0 


At  +  t, 


(div  V)2-  -^Ldiv  Vjt 


A  tJ  +  .  .  .  • 


} 


(26) 


Since  -jr-  div  V  =  0,  within  the  context  of  the  balance  equation  and  of  the 

order  (div  V)1  ,  the  quadratic  ttrm  of  eq.  (26)  approaches  the  linear  term  in 
moderate  latitudes  and  under  average  conditions  in  about  7  days.  The  replacement  of 
the  balance  equation  as  an  initialization  method  which  requires  a  moderate 
reprogramming  effort  becomes  thus  imperative  in  multi-level  numerical  weather 
prediction. 


According  to  Miyakoda  ct  al,37  the  movement  of  cyclones  and  the  tendency  for 
deepening  or  filling  are  the  major  problems  in  short-range  forecasts,  i.e.,  for  I  to  2 
days.  In  a  2- week  forecast,  the  life  histories  of  rycloncs  are  also  important  features  of 
the  prediction.  This  is  particularly  reflected  in  eq.  (26)  although,  for  medium-  to 
long-range  forecasts,  there  is  a  partial  adjustment  for  the  lack  of  an  initial 

div  V  =  0  because  of  the  interaction  of  the  primitive  equations,  the  continuity, 


and  the  thermodynamic  equation.  It  is  significant  and  in  correspondence  with  the 
above-slated  conclusions  that  in  the  experimental  predictions  described  by  Miyakoda 
et  al  the  intensities  of  the  highs  and  lows  weakened  appreciably  after  6  or  B  days 
reflecting  the  fact  that  the  forecast  of  eddy  kinetic  energy  was  less  than  the  observed. 
The  wiggling  in  the  pattern  of  geopotentiol  height  becomes  more  pronounced  with 
increasing  compulation  lime,  and  the  lack  of  development  of  a  certain  cyclone  on  the 


34  J.  Sma«orln*ky:  “On  the  Inclusion  of  Mol*t  Adiabatic  Procemc*  in  Numerical  Prediction  Model*." 

37  K.  Miyakoda  ct  tl:  “Experimental  Extended  Prediction*  with  a  Niiie-I.evel  llrmi*phrrlr  Model,”  Monthly 
Wwthnr  Rmilitw,  Vol  97,  No.  I,  Jan  I9ft9. 
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2d  and  3d  days  along  the  middle  Pacific  polar  frontal  zone  may  also  be  due  to 
inadequacies  of  initial  data  or  initialization  by  the  balance  and  associated  balance 
<u -equation.  The  quasi-stationary  modes,  or  long  waves,  are  more  dominant;  while  the 
eastward-moving  components,  the  relatively  shorter  waves,  are  too  small  in  amplitude 
which  is  in  agreement  with  Elsaesser’s  findings.  Significantly,  the  vertical  velocities 
calculated  by  the  balance  w  -equation  are  weaker  than  those  taken  from  the  prediction 
computation  based  on  the  time-dependent  primitive  equations. 

In  Sections  6  and  7,  it  becomes  apparent  that  the  filter  conditions  expressed  by 
eq.  (5)  are  cc  pletely  adequate  in  the  hydrostatic  system,  i.e.,  that  more  refined 
smoothing  techniques  utilizing  measured  winds  are  not  required  or  are  inconsistent 
with  the  hydrostatic  system  respectively.  Employing  the  hydrostatic  continuity 
equation  amounts  indeed  implicitly  io  a  smoothing  process  and  a  restriction  of  scale 
several  times  the  height  of  a  homogeneous  atmosphere  so  that  the  magnitude  of  the 
horizontal  divergence  does  not  exceed  that  of  the  vertical  vorticity.  Diffusion  terms, 
such  as  described  by  Smagorinsky  et  al,38  are  thus  not  due  to  a  lack  of  resolution 
inherent  in  the  filter  equations.  They  are  rather,  as  stated  by  the  authors,  a 
consequence  of  the  grid  size  used  which,  if  it  exceeds  30  km  significantly,  amounts  to 
an  additional  smoothing  'f  “signal”  functic  ts  obtainable  by  our  filter  equations. 

Finally,  it  should  be  mentioned  that  the  inequalities  (19)  and  (21 ),  which  are 
consistent  with  the  existence  of  the  continuity  equation  (10),  are  a  generalization  cf 
and  in  basic  agreement  with  the  results  obtained  by  Van  Mieghem,39  that  our  filter 
equations  are  sufficient  as  far  as  the  process  of  adaptation  described  by  Yeh  Tu-Cheng 
and  Li  Mai-Tsun40  is  concerned,  and  that  the  sigma-coordinate  system  developed  by 
Phillips41  is  also  subject  to  the  restrictions  of  the  hydrostatic  system  as  soon  as  the 
effective  static  stability  Hj  becomes  negative. 


f 

¥ 


6.  On  the  Modification  of  and  Initial  Conditions  for  the  Differential  Equations  of 
Meteorology  and  Related  Problems.  The  fundamental  problem  is  to  determine 
commensurable  smooth  fields  of  u,  v  and  w  under  consideration  of  the  fact  that  w 
cannot  be  measured,  generally  speaking.  Since  the  filtering  process  must  be  both 
consistent  for  all  three  wind  components  and  not  arbitrarily  independent  of  T.  p,  and 
p,  it  must  be  adaptive  in  nature.  Necessary  and  sufficient  conditions  for  an  optimal 
filtering  process  are,  therefore,  a  wind  vector  filter  operation,  complete  adapatability 
with  respect  to  the  continuity  and  thermodynamic  equation,  and  computability,  i.e., 
existence  of  a  unique  solution.  The  quasi-ergodic  filter  equation 


“d*  /  NblfLevei  Gener,d  Circulation  Model  of  the  Atmosphere,” 

Monthly  Weather  Review,  Vol  93,  No,  12,  December  1965. 


J.  M.  Van  Mieghem:  "Hydrodynamic  Instability,”  Compendium  of  Meteorology,  American  Meteorological 
Society,  Boston,  Massachusetta,  1951. 

40  Y.  Tu-Cheng  and  L  Mai-Taun:  "Adaptation  Between  the  Pressure  and  the  Wind  Field  in  Meao-  and  Small-Scale 
Motion*,”  AFCRL-68-0275,  Air  Forte  Cambridge  Research  Laboratories,  Bedford,  Massachusetts,  May  1968. 

N.  A.  Phillips:  "A  Coordinate  System  Having  Some  Special  Advantages  for  Numerical  Forecastim?/’  Journal  of 
Meteorology,  Vol  14,  Oct.,  1956. 
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which  yields  three  filter  conditions  involving  smooth  variables  u,  v,  and  &  fulfills  the 
requirements  mentioned  above  which  is  shorn  in  detail  in  Section  7.  It  has  to  be  kept 
d$ 

in  mind  that  is  to  be  expected  to  be  conserved  for  an  infinitesimal  part  of  the 


trajectory  only.  In  other  words:  The  resulting  three  diagnostic  filter  equations  arc 
more  accurate  than  eq.  (27)  implies. 


Inspection  of  the  filter  conditions  and  diagnostic  filter  equations  reveals  that 
transitory  sound  and  shear  gravity  waves  are  effectively  smoothed  out.  Typical  values 
for  the  wavelength  of  pressure  waves  are  100  m  to  1500  m,  and  shear  gravity  wave 
lengths  are  of  the  order  1000  m.  The  existence  of  these  waves  accounts  for  the 
quasi-iso  ti  opic  spectra  of  atmospheric  turbulence.  According  to  Charles41 ,  turbulence 
due  to  shear  in  the  vicinity  of  the  ground  exhibits  frequencies  above  2'A  cycles/second 
and  is  termed  “mechanical  turbulence”  in  contrast  to  “convection  turbulence”  which 
is  due  to  hydrostatic  instability  and  typically  occurs  at  frequencies  of  about  'A 
cycle/minute.  Since  the  analysis  of  shear  gravity  or  Helmholtz  waves  presupposes  the 
existence  of  wind  discontinuities  which  are  filtered  out  through  the  use  of  pressure 
and  temperature  fields,  a  wavelength  of  5  km  may  be  considered  as  the  lower  limit 
pertaining  to  the  applicability  of  the  new  filter  equations.  In  view  of  the  fact  that 
Anderson43  selected  a  300  m  height  interval  because  it  was  found  to  be  the  minimum 
interval  over  which  meaningful  values  of  both  vertical  wind  shear  and  lapse  rate  could 
be  derived  from  rawinsonde  data,  minimum  grid  distances  for  numerical  forecasts 
would  amount  to  300  m  —  400  m  in  the  vertical  and  about  2000  in  in  the  horizontal. 
Isolated  pressure  and  temperature  discontinuities  incompatible  with  such  three- 
dimensional  grid  would  have  to  be  smoothed  out  initially.  Indeed  application  of  the 
second-order  filter  condition 


d2& 

~dF 


=  0 


(28) 


and,  therefore,  abandonment  of  the  hydrostatic  assumption  is  a  priori  to  be  expected 
to  reduce  the  approximate  minimum  hydrostatic  wave  length  of  40  km  (five  times  the 
height  of  a  homogeneous  atmosphere)  by  one  order  of  magnitude.  In  agreement  with 
Pai44  we  can  thus  conclude  that  the  new  filter  equations  include  essentially  the 
vorticity  and  entropy  modes  of  the  hydrodynamic  equations,  but  not  the  transitory 


4iB.  N.  Charles:  “A  Summary  of  Available  Information  on  the  Vertical  Motions  of  Air  in  the  Troposphere,”  Sandia 
Corporation  Woridng  Paper  TID4500,  Office  of  Technical  Services,  Department  of  Commerce,  Washington, 
D.  C.,July,  1959. 

^A.  D.  Anderson:  “Free  Air  Turbulence, ”Joum.  of  Meteorology,  Vol.  14,  No.  6, 1957. 

^Shih-I  Pai:  "Viscous  Flow  Theory,"  Vol.  II,  D.  Van  Noetrand  Company,  Princeton,  NJ.,  1957. 
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sound  and  associated  shear  gravity  modes.  In  fact,  in  the  new  filter  equations  the 
Coriolis  terms  are  not  to  be  neglected  which  makes  them  amenable  to  meso-scale 
phenomena.  Only  small  scale  phenomena  such  as  individual  clouds  and  tornadoes 
cannot  be  incorporated. 


There  are  two  possibilities  as  to  the  use  of  filtered  wind  components  (t,  v,  (v  and 
corresponding  smooth  pressure  and  temperature  fields:  Either  these  fields  are  used  as 
quasi-exact  initial  values  together  with  a  very  small  grid  size,  horizontally  and 
vertically,  or  they  would  be  considered  as  approximate  initial  values  in  conjunction 
with  a  greater  grid  size.  In  the  first  case,  the  available  input  information  would  be 
utilized  to  its  maximum  extent,  and  commensurable  smaller  scale  “noise”  would  be 
generated.  In  the  second  case,  scale-preserving  operators  would  be  required,  either 
external  ones  such  as  a  two-dimensional  one  mentioned  by  Shuman45,  or  internal  ones 
of  the  diffusion  type.  Whether  a  small-scale  integration  can  be  achieved,  depends,  of 
course,  also  on  the  timely  availability  of  boundary  values.  We  do  not  adhere  to 
equations  involving  the  familiar  Reynolds  stresses  such  as 


aff  3u  .  a  3u 

-rr  +  -rr  +  u  r— 

3t  at  3x 


+ 


a  3u'  . 
11  ax 


(29) 


since  the  bars  lose  their  meaning  in  the  absence  of  an  crgodic  theorem  which  has  been 
emphasized  by  Kampe  de  Feriet.46  The  statistical  operation  shown  above  leads,  of 
course,  not  to  a  new  prognostic  equation.  For  example,  the  time  average  in  gen¬ 

eral  does  not  vanish  sufficiently  for  small  intervals.  Stable  statistical  averages  cannot 
even  be  obtained  over  longer  time  intervals  since,  as  is  well-known,  the  meteorological 
generation  process  is  not  a  stationary  one.  In  fact,  only  a  “potential  vorticity”of  the 
form 


f  +  *2*Bf4 


with  P  as  a  measure  of  standard  static  stability  and  ipB  as  the  stream  function 
obtained  through  the  classical  balance  equation  would  be  sufficiently  conserved.  For 
this  reason,  it  would  be  more  profitable  in  linear  statistical  forecasting  to  employ  \p- 
instead  of  ^-statistics.  The  use  of  bars  is  evidently  only  strictly  possible  in  the  case  of 
homogeneous  and  isotropic  turbulence  whence  the  wind-pressure  correlation  function 
vanishes.  In  addition,  incompressibility  is  assumed.  Under  these  conditions,  the 
Karman-Howarth  equation 


45 F.  G.  Shuman:  "Numerical  Method:  in  Weather  Prediction:  II  Smoothing  and  Filtering,”  Monthly  Weathei  llcview, 
Vol.  85,  No.  11. 

^J.  Kampe  de  Feriet:  Atmospheric  Diffusion  and  Air  Pollution.  Academic  Press,  New  York,  Unidon,  1959 
(Discussion,  p.  439). 
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2(ut)3/2  ^h'  +  4-^  =  2t>u 2  (f"+  4f')  (30) 

involves  only  wind  correlations.  In  view  of  the  above  it  must  be  stated  that 
fundamental  statistical  assumptions  associated  with  the  derivation  of  Reynolds  stresses 
do  not  hold  for  the  differential  equations  of  meteorology.  The  stresses  indicate  only 
that  the  generation  model  is  not  invariant  under  a  smoothing  operation  and  are  only  of 
use  in  connection  with  structural  and  analogy  considerations.  As  soon  as  the  link 
between  the  smoothed  and  perturbation  fields  is  assumed  to  be  negligible,  i.e.,  the 
generation  of  subgrid  noise  is  to  be  prevented,  diffusion  terms  have  to  be  added  or 
filters  have  to  be  applied  for  specified  time  intervals. 

In  view  of  an  adaptive,  i.e.,  coherent  smoothing  process  for  all  variables,  the 
horizontal  exchange  coefficient,  or  rather  function  K,  should  be  the  same  for  d,  0,  and 
fl.  The  diagnostic  filter  equation  for  (v  eliminates  the  need  of  diffusion  terms  for  this 
variable  and  identically  fulfills  the  equation  of  continuity  so  that  diffusion  of  mass 
does  not  occur.  According  to  Obuchow41  this  leads  to 

K  =  kiDo’3$a/3  (31> 

with  kj  as  a  dimensionless  factor  of  the  order  of  one.  23  as  an  observation  measure, 
and  Dq  as  proportional  to  the  squared  deformation  of  the  velocity  field 


Equation  (32)  is  only  an  approximation  since  it  has  been  derived  under  the  assumption 
of  homogeneous  turbulence  and  incompressibility,  whence  div  V  =  0.  It  holds, 
however,  well  in  its  two-dimensional  form  in  the  hydrostatic  systems,  whence  K 
=  k2  ^  .  In  this  form  it  has  been  utilized  by  Smagorinsky  et  a!.4®  The 

24/3  -relation  first  established  by  Richardson  has,  therefore,  to  be  considered 
questionable  for  regions  of  convective  activity.  This  has  been  ascertained  by  Lettau49 
and  is  due  to  strong  concurring  vertical  divergence.  Again,  only  because  of  the 
availability  of  the  diagnostic  and  simultaneously  optimal  (^-filter  equation  (63)  of 
Section  7,  eqs  (31)  and  (32)  in  their  two-dimensional  form  can  be  expected  to  hold 

4-A.  M.  Obuchow:  “Ueber  die  Energievcrtcilung  im  Spek'-'m  ciner  turbulr.ntcn  Stroemung,”  1941,  published  in 
Statistiche  Theorie  der  Turbulent,  Akademie-Veriag,  Berlin  1958. 

‘®J.  Smagorinsky,  et  al:  “Nine-Level  General  Circulation  Model  of  the  Atmosphere.” 

49H.  Lettau:  Atmospheric  Diffusion  and  Air  Pollution,  Academic  Press,  New  York,  London,  1959  (Discussion,  p. 
116). 
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reasonably.  In  case,  of  grid-scale  diffusion  involving  frequencies  which  are  small  in 
comparison  with  the  upper  limit  of  the  spectrum,  I)  *  const.  According  to 
Obuchow50  and  Richardson  eq.  (31)  reduces  then  to  k2  f(2).  This  simple  expression 
appears  to  apply  satisfactorily  to  geostrophic  conditions  only.  A.s  to  the  formulate  n 
of  suitable  vertical  diffusion  terms,  no  attempt  will  be  made  here. 


The  effects  of  ground  friction  are  not  considered  in  the  derivation  of  the  filtering 
equations  carried  out  in  Section  7.  Introduction  of  surface  stresses  of  the  form 


tx  ^pu/777  ,  ry  =  Cjjpv  Ju7 


(33) 
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lead,  however,  to  endogeneous  nonlinearities  with  associated  non-convergence  of  the 
numerical  relaxation  process.  Apparently,  quadratic  friction  terms  include  small-scale 
phenomena,  such  as  external  sound  and  shear  gravity  waves,  and  arc  also  not  strictly 
compatible  with  a  vanishing  vertical  velocity  at  a  horizontal  surface.  In  fact,  since  the 
variable  velocity  vectors  $(t)  at  the  upper  and  lower  boundaries  ar 1  not  available  and 
the  vertical  velocity  at  the  lower  boundary  is  a  function  of  an  approximate  horizontal 
velocity  dose  to  the  earth’s  surface  and  the  gcopotential  of  the  ground,  scales  below 
about  6  km  have  to  be  excluded  a  priori.  With  a  linear  frictional  term  introduced  by 
Guldberg  and  Mohn  and  proposed  as  suitable  by  Phillips51  the  first  equation  of  motion 
would  read  as 


du 

dt 


1 

'P 


dx 


+  fv  -  ku 


(3<i) 


|  More  explicitly,  k  =  a(x,  y)  k2(z)  though  k2  might  be  considered  constant  in  a 

{  sufficiently  deep  layer  adjacent  to  the  ground.  Application  of  terms  (33)  in  connection 

;  with  a  grid  of  small  length  would  only  lead  to  the  generation  of  undesirable  noise. 

1  Terms  of  the  form  ku  and  kv  can  be  fully  absorbed  in  the  filtering  equations  derived  in 

5  Section  7  although  they  have  been  omitted  for  the  purpose  of  simplicity.  Sin^e  motion 

in  the  friction  layer  is  such  that  there  is  horizontal  convergence  of  mass  in  areas  of 
j  cyclonic  activity  and  vice  versa,  as  stressed  by  Eliassen  and  Kieinschmidt,  these  terms 

j  ought  to  be  included  in  practical  computations. 

J 

j  The  incorporation  of  suitable  lower  and  upper  boundaries  in  our  initialization 

j  scheme  is  discussed  in  Section  7.  Pertaining  to  the  v/hole  globe,  the  problem  of  lateral 

t  boundary  conditions  does,  of  course,  not  present  itself.  In  a  quasi-hcmisphcric  grid 

system,  the  horizontal  wind  components  have  to  be  assumed  to  be  geostropldc,  and  an 

®A.  M.  OHuchow  “Ueber  die  Energieverlcilung  im  Spcktrum  eincr  turbulenten  Stroemung.” 

S,N.  A.  Phillip#:  Numerical  Weather  Prediction. 

52 A.  Eliaiwcn  and  E.  Kieinschmidt:  "Dynamic  Meteorology.” 
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approximate  vertical  velocity  has  to  be  computed  along  the  lines  described  in  Section 
5.  Lateral  winds  at  sizable  mountainous  shapes  have  in  general  to  he  assumed  as 
geostrophic  winds  also,  at  least  initially. 

The  Coriolis  parameter  has  been  treated  as  f  =  f(x,y)  in  Section  7.  For  the 
purpose  of  simplicity  we  have,  however,  considered  g  =  const.  In  a  more  exact  form, 
eq.  (43)  in  Section  7  would  read  as 


_d  3p  dp  dp 

dt  Sz  p  3F  « 


(35) 


“6> 

=  -~irw.  This  linear  correction  as  well  as  the  nonlinear  terms 


associated  with  and  ~  may  be  initially  neglected  in  the  relaxation  process, 


i.e.,  prior  to  obtaining  stable  numerical  results  from  the  simplified  as  well  as  linear 
filtering  equations. 


It  is  possible  to  derive  the  filtering  equations  for  spherical  coordinates  or  under 
consideration  of  map  scale  factors.  For  the  purpose  of  simplicity,  we  have  also  omitted 
diffusion  terms  in  the  continuity  equation  for  water  vapor  and  terms  involving 
differences  between  atmospheric  and  water  or  snow  temperatures  in  the  thermo¬ 
dynamic  equation.  These  complications  would  not  change  the  basic  character  of  the 
filtering  equations. 


The  discussion  in  the  preceding  paragraphs  appears  to  be  relevant  for  an 
understanding  of  the  ramifications  of  the  new  filtering  equations.  W?  have  not  only 
succeeded  in  deriving  optimal  filter  equations,  but  also  constructed  a  powerful  system 
of  four  prognostic  equations  and  one  diagnostic  equation  for  numerical  weather 
prediction  together  with  initial  conditions  the  details  of  which  are  exhibited  in  Section 
7.  This  system  would  sufficiently  cope  with  fronts,  hurricanes,  and  instabilities 
associated  with  the  jet  stream.  Were  it  not  for  the  inaccuracy  and  lack  of  observations 
and  the  restriction  of  an  upper  boundary  with  w  =  0,  this  system,  together  with 
well-digestible  non-adiabatic  inputs,  would  conceivably  make  it  possible  to  extend  the 
limit  of  predictability  associated  with  the  (x,  y,  p,  t)-system  considerably  beyond  2 
weeks.  This  may  be  comforting  to  scientists  involved  in  extended  forecasting  and 
general  circulation  simulations  such  as  Miyakoda53  and  Smagorinsky54  who  considered 
a  potential  limitation  of  2  weeks  as  rather  pessimistic. 


55  K.  Miyakoda,  et  il. 

u 

J.  Sm«£orinsky,  et  *1:  “Nine-Level  General  Circulation  Model  of  the  Atmosphere.” 


7.  Generalized  Optimal  Filter  Equations  Free  of  Hydrostatic  Limitations.  In 
Section  6,  we  have  shown  that  it  is  necessary  to  consider  pressure,  p,  a  priori  as  a 
continuous  variable  and  that  pressure  kinks  have  to  be  smoothed  out  in  order  to  avoid 
quasi-infinite  pressure  gradients  in  agreement  with  Haltiner  and  Martin.5*  Discontinu¬ 
ities  of  zeroth  order  involving  temperature,  or  rather  virtual  temperature,  require  even 
stronger  smoothing.  The  application  of  the  differential  equations  of  meteorology  is 
only  possible  with  smoothed  variables  including  consistently  filtered  winds.  For 
simplicity,  we  omit  the  filter  symbol  A  in  the  following  derivations. 

Since  the  filter  condition 

■ftr(pV)  =  0  (36) 

implicitly  includes  a  term  div  V ,  a  system  of  non-linear  partial  differential 

equations  would  result  which  does  not  permit  ar.  equilibrium  solution  and, 
consequently,  could  not  be  solved  by  relaxation  methods. 


As  already  mentioned  (Section  6),  the  filter  condition  is 


d*  V 
dtJ 


=  0 


(37) 


which  has  already  been  applied  in  Section  4  except  for  the  vertical  wind  component.  If 
eq.  (36)  were  valid,  it  would  imply  a  strong  identification  of  p  with  reference 
to  V  and,  thus,  make  the  differential  equations  more  deterministic,  i.e.,  all  invariant 
under  a  filter  operation,  which  is  not  possible. 


We  now  apply  eq.  (37)  with  reference  to  the  equations  of  motion 


du  X  aP 
dt  ”  p  3x 


+  fv 


(38) 


dv  1  3p 
dt  "p  3 y 


fu 


(39) 


dw  _  i  9p 
dt  ~  p  dz  ~  ® 


(40) 


55  G.  J.  Haltiner  and  F.  L  Martin:  Dynamical  and  Phytkal  Meteorology,  McGraw-Hill  Book  Company,  New  York, 
Toronto,  London,  19S7. 
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with  the  intermediate  result 

J  3p  1  dp  ,3p  d  , 

dt  ix  '  7  dt  '  $f  "P  dt  <fv>  =  0  (41) 

irly  “Tit-  *  (fu>  =  0  (42) 

d_9p  1  dp  3p  _  A 

dt  '  p  3t  *  a*  "  0  *  (43) 

With  F  =  i}£  under  consideration  of  eq.  (38)  and  (39)  with  regard 

*°  ~dt(^v)  *n  eq*  (41)  and  -jp  (f  v)  in  eq.  (42),  respectively,  and  in  view  of  the 

continuity  equation 

j  $  =  -  div  V  (44) 


we  arrive  at 


-H  |£  div  v  ~  (|H  |p  +  dw 

3x  \9x  3x  3x  3y  3x  )z ) 
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(  3f 
\  3x 
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(46) 

3F 

3z 

1 1£-  div  V  - 
tz 

/  3u 
\3z 

+ 

9z  3y  + 

toL  *2)  =  o 

3z  3z  j  u 

(47) 

In  the  not  step,  we  have  to  express  F  as  a  time-independent  function  which 
linearly  involves  the  divergence  f£  +  §}  +  The  thermodynamic  equation  reads 
in  a  very  general  form 

d(vf)  +  <«> 


*i**,*l®*'JI' “ ■  - ' •■  . ***“■ *  . — . . . *  —  . .  1  "  . . 


where  r,  designates  the  saturation  mixing  ratio,  L  *  600  cal  g' 1 ,  cp  =  .2405  cal  g" 1 
*  (deg C)*1  the  specific  heat  of  dry  air  at  constant  pressure,  c  =  1.0  cal  g'* ,  the  specific 
heat  of  water,  k  =  —  -  .2848  where  is  the  specific  heat  of  dry  air  for  constant 

volume,  pL  the  density  of  dry  air,  and  chj  non-precipitative  heat  added  to  a  unit 
mass  of  air.  With  a  =  -  ,  eq.  (48)  may  also  be  written  in  its  time-dependent  form 

cp 1 

ai^  +  (1.k.a  +  4,ri)d M-k^i  &. 

dlnr, 

Elimination  of  the  term  -  in  eq.  (49)  by  means  of  Smagorinsky’s  and  Col¬ 
lin’s  relation” 

dlnr,  dlnT  dlnp, 

It  Mt-D-jt-V1  <50> 

with  7  =  i  ftQtj'AftT  involves  L  as  the  latent  heat  of  condensation  and 

A'1  as  the  mechanical  equivalent  of  heat  leads  to 

[1  -  k  +u  .'7-2)  +  4.2 r,|  <k  +  a>  (51) 

From  eq.  (50)  and  (51)  follows 

dlnr,_  7k  +  a-l-4.2r$  din  pL 
dt  1  —  k  +  0(7  -  2)  +  4.2  r,  dt 

+ _ r _ Izl _ ^  .  J_  *1*  A  dln^  A  dlnp  52 

l-k  +  a(7-2)+4.2i;  c^T  dt  ¥  dt  ¥^T  {0^} 

Under  consideration  of 

V  [('  +  0.6r)  %■+  (1  +  0.6r) ,  O.fpT-^-j  (53) 

which  follows  from  the  equation  of  state,  in  view  of 


*  J.  Smaforirnky  and  G.  0.  Co  Ulna:  “On  the  Numerical  Prediction  of  Precipitation,”  Monthly  Weather  Review  83, 
1953. 
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with 


(55) 


{0  if  div  V  <  0  or  r  <  r, 

1  if  div  V  >  0  and  r  =  r, 
and  because  of 

ii=4i=0.6dr  <S6> 

Pl  P 

eq.  (53)  can  be  formulated  as 

?  =  (1+0-6'+4-S  •')  Tdi,V+fe7Tl' 

=  -  M[r,6]TdivV  +  N$L 


(57) 


Substitution  of  eq.  (57)  ineqs.  (45)  through  (47)  results  in  the  linear  diagnostic 
filter  equations 


MT 


feu+  +  [  mn. .  iLl  fei  +  ^  +is^ 

\3xJ  3x3y  bxbz)  [  3x  3xJ  \3x  3y  bz) 

A  au  JP  av  ap  aw  9p  f  A  fP  +  fu\ 

3x  3x  3x  5^  3x  3z  ”  \p  by  ) 

♦4£  -If) 


(58) 


MT/  32u  ,  32 v  ■  9aw\  f.  [ 3(MT)  3p]  Au  +  *1  + 

MT\3x3y  +  Sy*  bybz)  [  9y  <^y J  9y  ^iz) 

,3u3T>+iv9p  3w9p  +fflf_L£P  -  fv\ 

3y  3x  3y  3v  3y  3z  p\p  bx  J 

.Up(u|I  +  v|I)  (n^)  =  o 


(59) 
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We  have  to  remember  that  the  saturation  mixing  ratio  r,  =  rt(p,  T)  and  7  =  7(T) 
and  that  it  is  necessary  to  obtain  first  a  solution  of  eqs.  (58),  (59)  and  (63)  with  8  =  0 
whereupon  the  criterion  (55)  is  applied.  One  or  two  iterations  will  then  yield  satisfac¬ 
tory  results.  Unless  6  has  some  variability  along  the  lines  suggested  by  Smagorinsky,57 
the  variable  5  should  be  about  0.8  instead  of  1.0  in  agreement  with  numerical 
simulations. 

It  is  to  be  expected  that  the  over-relaxation  factors  d  in  the  iteration  scheme 

U(n  +  D  =  u$»  +  a  G,  [u(n)  v(n),  w(n)j 

v<n+1>  =  v<">  +  0  Gj  [u(n),v(n\  w(n>]  (66) 

with  u<‘>  =  -  ^  ,  V<1>  =  £  .  *(1>  =  0 

in  which  Gj  and  G2  represent  residuals  of  eq.  (58)  and  eq.  (59),  respectively,  have 
to  be  quite  small  in  small-scale  solutions  involving  strong  divergence  (convergence)  and 
vertical  wind  velocities.  Since  u*!*  and  v(1)  become  singular  at  the  equator  and 
convergence  is  slow  in  very  low  latitudes,  fine  grid  solutions  are  not  possible  in  the 
vicinity  of  the  equator.  Due  to  the  fact  that  the  mass  field  cannot  be  accurately 
determined  in  the  equatorial  region,  horizontal  winds,  obtained  through  the  tracking 
of  floating  balloons,  and  additional  temperature  measurements  would  facilitate  the 
computation  of  all  desired  quantities.  The  use  of  diagnostic  filtering  equations  for  this 
purpose  has  been  mentioned  by  several  authors  including  Mintzsi  though  in 
connection  with  the  more  restrictive  hydrostatic  prediction  system. 

Utilization  of  the  hydrostatic  approximation  with  respect  to  height  determina¬ 
tions  weakens  the  application  of  the  filtering  and  associated  prognostic  equations  as 
far  as  smaller  scales  are  concerned  but  still  allows  the  computation  of  divergences 
exceeding  the  vorticity  on  a  constant  pressure  surface.  This  is  of  importance  pertaining 
to  the  immediate  applicability  of  the  new  prediction  system. 

As  to  the  upper  boundary  condition,  the  assumption  w  =  0,  of  course,  has  to 
be  made  for  a  finite  height.  In  this  inspect,  the  condition  var  w  =  Min.  would  provide  a 

good  separation  criterion.  This  has  to  coincide  with  the  criteria  var  -|^-  =  Min., 

var  =  Min.,  var  =  Min.,  var  =  Min.,  as  far  as  interpolation  from  a 

lower  to  a  higher  level  is  concerned,  i.c.,  to  an  average  equilibrium  boundary  which 

nY.  Mintz;  “The  Four  Baric  Requirement*  for  Numerical  Weather  Prediction  in  Global  Weather  Prediction,"  ed.  by 
Bruce  and  KieJy,  Holt,  Rinehart  and  Winston,  New  York,  1970. 

!*J.  Smagorinsky:  "On  the  Dynamical  Prediction  of  Large-Scale  Condensation  by  Numerical  Methods,” 
Geophysical  Monograph  No.  S,  American  Geophysical  Union,  1960. 
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exists  at  about  20-km  height.  In  addition,  linear  correlation  based  on  considerable 
prior  work  would  permit  the  extrapolation  of  u,  v,  p,  and  T  and  thus  facilitate  a 
w-caleulation  in  virtue  of  eq.  (63). 


U.  Concluding  Remarks.  It  is  important  to  emphasize  that  initialization  and  tin; 
structure  of  prognostic  equations  are  intimately  related.  Furthermore,  the  determina¬ 
tion  of  the  w-field  must  be  consistent  with  that  of  the  u-  and  v-fields.  If  the 
differential  equations  of  meteorology  were  to  be  applied  in  a  eompletely  unmodified 
form,  including  all  equations  of  motion,  the  initial  fields  would  have  to  be  found 
simultaneously  in  a  grid  with  very  small  lengths  Ax,  Ay,  Az  which  is  impossible.  The 
optimal  filter  equations  are  also  conditional  equations  with  reference  to  the  wind 
fields  once,  sufficiently  smooth  pressure,  temperature  and  humidity  fields  are  given. 
Since  these  equations  imply  a  very  delicate  “balance,”  independently  given  fields  of  all 
variables  excluding  the  generally  not  available  w-field  would  lead  to  the  generation  of 
artificially  large  gravity  oscillations  as  stressed  by  Phillips.59  Because  of  the  existence 
of  thermal  convection,  the  hydrostatic  differential  equations  cannot  adequately 
describe  the  physical  processes  in  sufficient  detail  and  appear  to  be  unsuitable  as  a 
basis  for  parameterization  processes  including  cumulus  convection  which  have  been 
emphasized  by  Leith.*0  In  view  of  the  above-mentioned  facts,  the  non-availability  of 
timely  boundary  values  and  stability  problems,  a  practical  forecast  system  of  highest 
predictability  is,  consequently,  one  which  contains  only  one  diagnostic  equation.  The 
existence  of  such  an  equation  is  equivalent  with  the  existence  of  consistent  initial 
fields  of  variables.  The  initialization  and  prognostic  model  developed  here  has  a 
considerably  greater  predictive  ability  and  versatility  than  the  equations  presently  in 
use.  Long-range,  nun-adiabatic  effects  other  than  precipitation  would  also  be  of  greater 
significance  in  context  with  the  new  system  of  equations.  These  equations  arc  also 
indispensable  with  respect  to  the  dynamical  prediction  of  medium  to  large-seale 
condensation  and  many  research  efforts  including  general  circulation  simulations.  As 
can  be  ascertained  from  a  number  of  statements  made  by  the  Committee  on 
Atmospheric  Sciences,  National  Research  Council,  and  listed  in  the  appendix,  this 
research  has  made  contributions  to  a  variety  of  outstanding  problems.  Finally,  this 
study  also  has  implications  regarding  kind,  methods,  density,  and  frequency  of 
measurements  and  indicates  the  necessity  for  a  greater  effort  of  and  cooperation 
between  the  National  Oceanic  and  Atmospheric  Administration  and  the  Aeronautics 
and  Space  Administration  and  beyond  that,  internationally. 


59N.  A.  Phillips:  “Numerical  Weather  Prediction." 

*°C.  Itflth:  “A  Six  Level  Model  of  (lie  Almoaphere"  In  "Global  Weather  Prediction,"  rd.  by  llmee  and  Klely,  Holt, 
Ilyich  art  v‘>d  Wlnxlon,  New  York,  1970 
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APPENDIX 


RELEVANT  STATEMENTS  OF  THE  COMMITTEE  ON 
ATMOSPHERIC  SCIENCES,  NATIONAL  RESEARCH  COUNCIL 

The  following  statements  on  outstanding  problems  in  weather  prediction  pertain¬ 
ing  to  which  contributions  have  been  made  in  this  research  have  been  published  in 
“The  Atmospheric  Sciences  and  Man's  Needs,”  National  Academy  of  Sciences,  Wash¬ 
ington,  D.  C.,  1971. 

1.  Synoptic  and  Planetary  Scales,  p.  18,  19.  ...  Progress  has  been  limited  due  to 
the  computational  requirements  and  to  the  fact  that  vertical  convection  and  condensa¬ 
tion  play  such  important  parts  in  mesoscale  phenomena.  Although  numerical  models 
have  successfully  simulated  convection  cells,  the  dynamics  and  statistics  of  convection 
have  not  been  incorporated  in  the  models  of  the  large-scale  circulation;  this  constitutes 
one  of  the  most  serious  limitations  on  progress  in  atmospheric  prediction.  . . . 

2.  Changes  of  Microscale  and  Mesoscale,  p.  19-21.  . . .  There  are  severe  mathemati¬ 
cal  and  theoretical  difficulties  in  developing  general  prediction  models  for  these  smaller 
scales.  The  time  scales  are,  of  course,  so  short  that  the  quasi-geostrophic  relation  be¬ 
tween  wind  velocity  and  pressure  is  not  valid.  Furthermore,  many  of  the  weather  sys¬ 
tems  are  fully  three  dimensional,  so  that  the  hydrostatic  approximation  which  is  essen¬ 
tial  to  general  circulation  theory  cannot  be  applied  to  the  smaller  scales.  . . . 

. . .  Other  phenomena  of  the  mesoscale  and  microscale,  e.g.,  vertical  convection, 
clear-air  turbulence,  and  secondary  features  of  the  planetary  boundary  layer,  have  not 
been  incorporated  into  the  numerical  models.  The  physics  of  these  weather  features, 
especially  their  interaction  with  phenomena  of  other  scales,  remains  largely  unknown, 
and  promising  ideas  are  needed.  . . . 

3.  Recommendations,  p.  23-28.  .  .  .  Scientific  plans  for  the  U.  S.  contribution  to 
GARP  were  outlined  in  the  report  of  the  U.  S.  Committee  for  GARP  referred  to  earlier. 
Bhat  report  emphasized  the  following  requirements:  .  .  .  (c)  the  conduct  of  regional 
field  programs  and  computer  modeling  experiments  to  improve  the  physical  and  mathe¬ 
matical  basis  of  long-range  prediction.  . . . 

.  .  .  There  are,  however,  several  specific  problems  that  have  been  identified  by 
USC-GARP  and  that  require  prompt  attention  in  order  that  the  GARP  can  be  success¬ 
fully  carried  out.  These  are  .  .  .  (b)  A  strong  numerical  modeling  effort  is  needed  to 


carry  out  simulation  experiments  aimed  at  design  of  the  global  observing  system.  . . . 

. . .  Improvements  in  prediction  models  and  the  associated  increased  understand¬ 
ing  of  atmospheric  processes  will  also  contribute  to  the  necessary  basis  for  virtually 
every  other  application  considered  in  this  report.  . . . 

. . .  Recommendation  1-4.  In  order  to  advance  understanding  of  mesoscale  phe¬ 
nomena  and  to  improve  ability  to  forecast  these  phenomena,  emphasis  should  be 
placed  on  research  on  fronts,  jet  streams,  organized  convection,  and  other  phenomena 
of  mesoscale. 

Increased  support  should  be  given  to  the  development  of  numerical  models 
capable  of  predicting  such  features  as  the  weather  associated  with  fronts  and  oiganized 
convection.  The  largest  available  computer  in  the  hands  of  a  competent  numerical 
modeling  group  is  needed  for  this  purpose.  . . . 
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